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ABSTRACT 

Arguments for and against the widely accepted picture of a solar dynamo being seated in the 
tachocline are reviewed and alternative ideas concerning dynamos operating in the bulk of the con- 
vection zone, or perhaps even in the near-surface shear layer, are discussed. Based on the angular 
velocities of magnetic tracers it is argued that the observations are compatible with a distributed 
dynamo that may be strongly shaped by the near-surface shear layer. Direct simulations of dynamo 
action in a slab with turbulence and shear are presented to discuss filling factor and tilt angles of 
bipolar regions in such a model. 

Subject headings: MHD - Sun: magnetic fields - sunspots 



1. INTRODUCTION 

There appears to be general consensus that the solar 
magnetic field is generated and stored in the overshoot 
layer near the bottom of the convection zone (Spiegel 
& Weiss 1980, Golub et al. 1981, Galloway & Weiss 
1981, van Ballegooijen 1982, Choudhuri 1990). This 
layer is now believed to coincide with the tachocline, i.e. 
the radial shear layer at the bottom of the convection 
zone where the latitudinal differential rotation changes 
into rigid rotation in the radiative zone (Spiegel & Zahn 
1992). The main arguments in favor of this proposal are 
connected with flux storage over times long enough for 
the shear to amplify the toroidal field (Moreno-Insertis, 
Schiissler, & Ferriz-Mas 1992), and with the observed 
size of active regions (~ 100 Mm) being comparable with 
the typical eddy scale at the bottom of the convection 
zone (Galloway & Weiss 1981), as well as the observed 
fidelity of Hale's polarity law; see Fisher et al. (2000), 
Tobias (2002), Schiissler (2002), Ossendrijver (2003), Fan 
(2004), and Weiss (2005) for recent reviews. All these as- 
pects are intimately related to the thin flux tube picture. 
Indeed, one of the big successes of the thin flux tube ap- 
proximation (Spruit 1981, Moreno-Insertis 1986, Chou & 
Fisher 1989) is the quantitative prediction of Joy's law 
describing the latitudinal dependence of the observed tilt 
angles of bipolar regions. It is found that Joy's law is 
obeyed only for flux tubes with magnetic fields that are of 
the order of 10 5 G (D'Silva & Choudhuri 1993, Schiissler 
et al. 1994, Caligari, Moreno-Insertis, & Schiissler 1995). 
This result poses rather stringent demands on dynamo 
theory that are hard to meet. Although it may already 
be hard for the differential rotation in the tachocline to 
amplify a poloidal field to a strength of ~ 10 5 G, which 
may require flux intensification by exploding flux tubes 
(Rempel & Schiissler 2001), it is not obvious how to ex- 
plain the production of strong and sufficiently coherent 
poloidal field that is necessary to produce the toroidal 
field. 

The purpose of the present paper is to discuss the diffi- 
culties for dynamo theory in meeting these demands and 
to reconsider the alternative scenario that the solar dy- 
namo may operate in the bulk of the convection zone, or 



perhaps even in the near-surface shear layer in the upper 
35 Mm of the sun. This is the layer where recent he- 
lioseismological inversions have shown marked negative 
radial shear (Howe et al. 2000a, Corbard & Thompson 

2002, Thompson et al. 2003). The presence of a deeper 
layer that spins about 5% faster than the photosphere 
has always been anticipated based on the higher rota- 
tion rate of magnetic tracers (Gilman & Foukal 1979, 
Golub et al. 1981). It remained unclear, however, just 
how deep or shallow this layer really is. The most nat- 
ural assumption at the time was to place this layer near 
the bottom of the convection zone where magnetic buoy- 
ancy is weak and shear could work on the field unim- 
peded by the turbulence. In the early days of mean 
field dynamo theory a negative radial f2 gradient was al- 
ready anticipated (Stix 1976, Yoshimura 1976) because 
it would explain the observed anticorrelation of the signs 
of mean azimuthal field (inferred from the orientation of 
bipolar regions) and the mean radial field (measured by 
magnetograms). This is because negative radial shear 
turns a positive radial field into a negative azimuthal 
field, producing the observed anticorrelation. A negative 
radial gradient of angular velocity seemed confirmed by 
observations of very young sunspots that rotate faster 
than older ones (Tuominen 1962), suggesting that they 
may be anchored in the layer where the angular veloc- 
ity is maximum (Tuominen & Virtanen 1988, Balthasar, 
Schiissler, & Wohl 1982, Nesme-Ribes, Ferreira, & Mein 
1993, Pulkkinen & Tuominen 1998). 

The sunspot observations are not easily explained by 
interface dynamos, unless one is able to show that the 
angular velocity of magnetic tracers is just the pattern 
speed of a traveling wave phenomenon. A similar pro- 
posal has been made to explain the observed pattern 
speed of the supergranulation (Gizon, Duvall, & Schou 

2003, Schou 2003, Busse 2004). 

The importance of the near-surface shear layer has al- 
ready been investigated by Dikpati et al. (2002) who 
studied the effects of near-surface radial shear on a flux 
transport dynamo. They came to the conclusion that the 
effect of the near-surface shear layer is subdominant in 
the context of the flux transport model studied earlier 



2 



by Dikpati & Charbonneau (1999). Mason, Hughes, & 
Tobias (2004) also considered the issue of near-surface 
dynamo action in the context of a two-layer dynamo 
(one at the top and one at the bottom of the convec- 
tion zone). They allowed for an additional a effect in 
the upper layers, retaining only the radial shear in the 
tachocline. They came to the conclusion that the near- 
surface dynamo was harder to excite due to the larger 
distance to the tachocline. 

There are several arguments in favor of a dynamo oper- 
ating in or being strongly controlled by the near-surface 
layer of the sun. First, in the outer 35 Mm (correspond- 
ing to fractional radii r/R > 0.95) the negative radial 
shear, together with an a effect of the usual sign, would 
easily explain the observed equatorward migration of 
sunspot activity (e.g. Parker 1979). Second, the nega- 
tive phase relation between radial and azimuthal mean 
fields, B r and B^, respectively, would be automatically 
satisfied (Stix 1976, Yoshimura 1976). Third, the ra- 
dial near-surface shear is particularly strong between the 
equator and 30° latitude, and weak near the poles; see 
Fig-E ln th° tachocline, by contrast, there is hardly any 
radial shear at 30° latitude and maximum shear near the 
poles. Invoking near-surface shear for producing sunspot 
activity would naturally explain the strongest production 
of sunspot activity at low latitudes and the much weaker 
magnetic activity at high latitudes with a possible pole- 
ward migration (Stix 1974), provided the near surface 
shear changes sign at high latitudes, as is perhaps indi- 
cated by helioseismology (Thompson et al. 2003). [The 
poleward branch could also be due poleward flows - as 
demonstrated convincingly by flux transport models; see 
Baumann et al. (2004) and Schussler (2005).] Fourth, 
the rotational velocity of very young sunspots (age less 
than 1.5 days) is 14.7°/day at low latitudes (Pulkkinen 
& Tuominen 1998), corresponding to 473 nHz, which is 
about the largest angular velocity measured with helio- 
seismology anywhere in the sun; cf. Fig. ^ This corre- 
sponds to the angular velocity at a radius r/R = 0.95, 
which is 35 Mm below the surface. Similar conclusions 
can be drawn from the apparent angular velocity of old 
and new magnetic flux at different latitudes (Bcncvolcn- 
skaya et al. 1999). We will return to these observations 
in 

The main argument against distributed and near- 
surface shear layer dynamos is that magnetic flux tubes 
are highly buoyant in the convection zone proper (Parker 
1975). Thus, too much magnetic flux may be lost, on 
time scales so short that the shear cannot amplify the 
poloidal field to significant field strengths. However, over 
the past 15 years simulations have shown that magnetic 
buoyancy is strongly offset by the action of turbulent 
pumping, which leads to a net accumulation of magnetic 
energy at the bottom of the convection zone (Nordlund 
et al. 1992, Brandenburg et al. 1996, Tobias et al. 1998, 
2001, Dorch & Nordlund 2001, Ossendrijver et al. 2002, 
Ziegler & Riidiger 2003). Indeed, magnetic pumping 
is now also being invoked to keep down the horizontal 
magnetic fields just outside the penumbra of sunspots 
(Thomas et al. 2002, Weiss et al. 2004). Thus, we can 
envisage pumping as a mechanism that tries to keep the 
surface clean of magnetic fields, but it can do so only 
approximately and only if the field is not already too 
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Fig. 1. — Radial profiles of the internal solar rotation rate, 
as inferred from helioseismology. The angular velocities of active 
zones at the beginning of the cycle (at ?s 30° latitude) and near 
the end (at fa 4°) is indicated by horizontal bars, which intersect 
the profiles of angular velocity at r/R a 0.97. [Adapted from 
Benevolenskaya et al. (1999).] 



strong. 

Another argument against near-surface dynamos is the 
high degree of turbulence in the upper layers, which could 
lead to strong random distortions of the orientation of 
flux tubes. This leads to the notion of rising flux tubes 
being 'brain washed' during their ascent (Schussler 1983, 
1984), i.e. they lose their original east-west orientation 
and would not obey Hale's polarity law. However, this 
picture derives originally from the idea that flux tubes 
are produced in deep layers at or below the overshoot 
layer and are then subjected to a more passive buoyant 
rise through the convection zone. Here, however, we are 
envisaging the production of sunspots much closer to the 
surface, so the notion of flux tubes rising through a ma- 
jor portion of the convection zone is not invoked. Indeed, 
local helioseismology suggests a picture quite compat- 
ible with sunspots being a shallow surface phenomenon 
(Kosovichev, Duvall, & Scherrer 2000, Kosovichev 2002). 
The actual sunspot formation might then be the result 
of convective collapse of magnetic fibrils (Zwaan 1978, 
Spruit & Zweibel 1979), possibly facilitated by negative 
turbulent magnetic pressure effects (Kleeorin, Mond, & 
Rogachevskii 1996) or by an instability (Kitchatinov & 
Mazur 2000) causing the vertical flux to concentrate into 
a tube. 

It should be noted that the picture of shallow sunspots 
does not necessarily contradict the idea of strong flux 
tubes rising to the surface. In fact, as the tube rises to the 
surface, it must eventually undergo catastrophic expan- 
sion (Moreno-Insertis, Caligari, & Schussler 1995). This 
would detach the forming active region and its sunspots 
from its roots (Schrijver & Title 1999, Schussler 2005), 
which might then be compatible with the shallow sunspot 
picture from local sunspot helioseismology. 

Yet another potential problem with near-surface shear 
layer dynamos are the relatively short turbulent time 
scales. However, in 35 Mm depth the typical turnover 
time is, according to mixing length theory (Spruit 1974), 
already 1-3 days. Therefore the inverse Rossby number, 
2£1t, is of the order of unity, so the turbulence is certainly 
beginning to be affected by rotation. As is familiar from 
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mean field dynamo theory, the combination of poloidal 
and toroidal fields really corresponds to a right-handed 
spiral in the northern hemisphere. Thus, whenever parts 
of this spiral touch the surface they produce a bipolar 
active region with the observed tilt angle. This will be 
discussed further in H4.4I 

We now turn to the question of small scale magnetic 
fields. For magnetic Prandtl numbers as small as those in 
the sun, the magnitude of the turbulent magnetic fields 
from local small scale dynamo action at the top of the 
convection zone (Cattaneo 1999) is possibly not much 
stronger than the magnitude of fields from the large scale 
dynamo. This suggestion is motivated by the recent re- 
alization that small scale dynamo action (as originally 
explored by Kazantsev 1968) becomes either completely 
impossible or at least much harder to excite when the 
magnetic Prandtl number becomes small (Schckochihin 
et al. 2004a, Boldyrev & Cattaneo 2004, Haugen, Bran- 
denburg, & Dobler 2004). The dominance of small scale 
over large scale dynamo activity in global simulations 
of solar-like convection (Brun, Miesch, & Toomre 2004) 
might therefore also be related to the fact that the mag- 
netic Prandtl number is not small in those simulations. 

Another important aspect is the fact that in the pres- 
ence of shear, turbulent dynamos can produce and main- 
tain fields of equipartition strength (Brandenburg et al. 
2005). We will discuss some of those models also in 
2] Although such models still lack important aspects 
of solar dynamos (convection, stratification, and rota- 
tion), they are quite suitable for testing new effects in 
mean field theory, for example the shear-current effect 
(Rogachevskii & Kleeorin 2003, 2004) and current he- 
licity fluxes (Brandenburg & Sandin 2004, Subramanian 
& Brandenburg 2004). We begin by discussing first the 
problems associated with the current picture of dynamos 
operating in the tachocline. A summary of arguments 
discussed in the text is given in Tabled 

2. PROBLEMS WITH TACHOCLINE DYNAMOS 

By tachocline dynamos we mean dynamos where the 
main shear that is responsible for the cyclic toroidal fields 
originates from the tachocline. These dynamos take into 
account the measured differential rotation profile, al- 
though sometimes the latitudinal shear is neglected (e.g. 
Parker 1993, Choudhuri, Schiissler, & Dikpati 1995). 
Tachocline dynamos can be divided into three main sub- 
classes: (i) overshoot dynamos where there is only a neg- 
ative a effect in the overshoot layer, (ii) interface dy- 
namos where a negative a effect is assumed in the upper 
parts of the convection zone, and (iii) Babcock-Leighton 
type flux transport dynamos where the a effect is also lo- 
cated near the surface, but it is now positive and there is 
meridional circulation transporting flux in the overshoot 
layer from high latitudes toward the equator. 

One of the longstanding problems with dynamos oper- 
ating in a thin layer at the bottom of the convection zone 
is the large number of oppositely oriented toroidal flux 
belts in each hemisphere (Moss, Tuominen, & Branden- 
burg 1990). This tends to produce a rather unrealistic 
butterfly diagram. This problem can partly be allevi- 
ated by increasing the thickness of the overshoot layer 
(Riidiger & Brandenburg 1995) to about 50 Mm, which 
is beyond the currently accepted thickness of the over- 
shoot layer of about 7 Mm or less (Basu 1997). 



Another problem is the strong radial shear at polar 
latitudes in the tachocline. This leads to a dominance of 
magnetic activity in polar regions. It is therefore custom- 
ary to postulate an artificially modified latitudinal de- 
pendence of the a effect. Riidiger & Brandenburg (1995) 
assumed that a was proportional to cos 2 #, where 6 is 
colatitude, and Markiel & Thomas (1999) assumed that 
a was proportional to a gaussian concentrated around 
the equator. This manipulation was originally motivated 
by the possible presence of higher order terms quantify- 
ing the combined influence of stratification and rotation. 
Simulations of Ossendrijver et al. (2002) have indeed con- 
firmed a suppression of a near the poles. Nevertheless, 
it remains puzzling that at 30° latitude, where sunspots 
first emerge, the radial shear in the tachocline basically 
vanishes. So there should not be any local toroidal field 
generation. This problem may however be alleviated in 
the context of flux transport dynamos, as will be dis- 
cussed at the end of this section. 

We recall that the positive radial angular velocity gra- 
dient in the tachocline stretches a positive B r into a pos- 
itive B^,, so that B r B,p is also positive. As discussed 
in the introduction, this is in conflict with observations 
(Yoshimura 1976, Stix 1976). Although in some models 
B r B^ can still be negative during certain intervals and 
in certain latitudes, this cannot be regarded as a robust 
or well understood feature. Also, the occasional intervals 
of positive B r B^ seen in some models (Kiiker, Riidiger, 
& Schultz 2001) depend on assumptions about the depth 
were the toroidal field is evaluated. Furthermore, these 
models rely on the negative a that is expected at the 
bottom of the convection zone (Yoshimura 1972, Krivo- 
dubskii 1984, Riidiger & Kitchatinov 1993). 

If the cyclic field in the sun really originates from the 
tachocline, one would expect to see cyclic modulations 
of the local angular velocity, similar to those seen at the 
surface (Howard & LaBonte 1980). Instead, there is pos- 
sible evidence for a shorter 1.3 yr period at the base of the 
convection zone (Howe et al. 2000a) . This would suggests 
that the field responsible for the 22 yr cycle cannot come 
from the tachocline, but rather from the outer 70 Mm 
of the sun where 11 yr variations have indeed been seen 
(Howe et al. 2000b, Vorontsov et al. 2002). Indeed, a 
recent model by Covas et al. (2001) explains the 1.3 yr 
period in terms of spatio-temporal fragmentation, where 
the dynamo has a shorter period at the bottom of the 
convection zone and a longer period in the upper parts. 
In this model the field responsible for the 22 yr cycle 
would originate from the upper parts of the convection 
zone (CZ). 

We note in passing that, if a tachocline was really cru- 
cial for a dynamo to work, one might expect a break in 
the magnetic activity toward late M dwarfs that become 
fully convective and therefore lack a tachocline. This 
is not observed (Vilhu 1984). However, this argument 
is not really compelling because it can be argued that 
the fields of fully convective stars are only of small scale 
(Durney, De Young, & Roxburgh 1993). 

Finally, a problem with any model drawing its field 
from deep underneath is that it is not easy to imagine 
that a flux tube can maintain its integrity while rising 
over 20 pressure scale heights from the bottom of the 
convection zone to the top. Indeed, direct simulations 
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TABLE 1 

Summary of arguments for and against tachocline and distributed dynamos discussed in the text 



arguments 



tachocline dynamos 



distributed/near-surface dynamos 



in favor 



against 



flux storage 

turbulent distortions weak 

correct butterfly diagram with meridional circulation 
size of active regions (~ 100 Mm) naturally explained 



100 kG field hard to explain 
flux tube integrity during ascent 
too many flux belts in latitude 
maximum radial shear at the poles 
no radial shear where sunspots emerge 
quadrupolar parity preferred 
wrong phase relation 

1.3 yr instead of llyr variation of at base of CZ 
coherent meridional circulation pattern required 



negative surface shear yields equatorward migration 
correct phase relation 

strong surface shear at latitudes where the spots are 

max(C)/27r = 473 nHz agrees with f2(youngest spots) 

active zones move with £7(0.95) 

llyr variation of f2 seen in the outer 70 Mm 

even fully convective stars have dynamos 

strong turbulent distortions 
rapid buoyant losses 

too many flux belts if dynamo only in shear layer 

not enough time for shear to act 

long term stability of active regions 

profile of Q(youngest) by 4 nHz above Q(0.95) 

possible anisotropies in supergranulation 



show that a large amount of twist is needed to keep the 
tubes intact over at least a few pressure scale heights 
(Moreno-Insertis & Emonet 1996). On the other hand, 
although a modest amount of twist can be useful for 
explaining the so-called S spots (Fan et al. 1999), too 
much twist can make the tubes kink- unstable (e.g., Lin- 
ton, Longcope, & Fisher 1996). It is also not clear how 
the faster sunspot proper motion of very young sunspots 
(Pulkkinen & Tuominen 1998) can be explained if the 
spots were rooted in the tachocline. It would be much 
more straightforward if they were rooted near the maxi- 
mum of fl about 35 Mm below the surface. 

A different class of models are the Babcock-Leighton 
type flux transport dynamos, as recently studied by Dik- 
pati & Charbonneau (1999), Nandy & Choudhuri (2002), 
and Chatterjee et al. (2004). Here the a effect is assumed 
to come from the surface layers, so a would be positive. 
These models deal with some of the aforementioned prob- 
lems by invoking a grand meridional circulation pattern 
playing the role of a conveyor belt that transports flux 
through the tachocline from high latitudes to the equa- 
tor. This circulation is responsible for driving the dy- 
namo waves equatorward and are also determining the 
cycle period (Durney 1995, Choudhuri et al. 1995, Kiiker 
et al. 2001). A polar branch, on the other hand, can be 
explained by postulating a two-cell circulation pattern. 
It is this type of model for which the effect of the near- 
surface shear layer has recently been investigated by Dik- 
pati et al. (2002). However, form and magnitude of the 
meridional circulation in the sun are quite uncertain. Di- 
rect simulations by Miesch et al. (2000) suggest a rather 
more irregular pattern of many cells changing with time. 
If this result continues to persist also in more realistic 
simulations, it would render the flux transport picture 
rather fragile. 

3. PROBLEMS WITH DISTRIBUTED DYNAMOS 

We begin by discussing first the evidence from mag- 
netic tracers in favor of their near-surface anchoring and 
turn then to the more theoretical arguments supporting 
the notion of distributed dynamos that are being strongly 
affected by the near-surface shear layer. 
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Fig. 2. — Angular velocity of sunspots of different age as a func- 
tion of colatitude 8 compared with the helioseismological internal 
angular velocity at r = 0.95-R and with the Doppler velocity at the 
surface. The large gray dots at about 30° and 4° denote the an- 
gular velocities of active zones respectively near the beginning and 
near the end of the cycle. The three small dots near the zig-zag 
line are the helioseismological angular velocities at r = 0.7R. 



Magnetic tracers have long been known to rotate faster 
than the photospheric plasma (e.g. Gilman & Foukal 
1979). One possibility is that magnetic flux tubes pos- 
sess an angular velocity that is in excess of that of the 
surrounding plasma (Wilson 1987). Using the data of 
Pulkkinen & Tuominen (1998) we plot in Fig. [3 the an- 
gular velocity of sunspots of different age (from 1.5 days 
to several months) versus colatitude, where we have fit- 
ted their angular velocity to the common representation 
fl = a + bcos 2 8, which leads to a linear graph when 
plotted versus cos 2 6. These results are compared with 
the helioseismologically determined angular velocity at 
r = 0.95i? and 0.7i?, as well as with the Doppler veloc- 
ity at the surface, using a fit up to cos 4 9, as quoted by 
Thompson et al. (2003). 

Let us now discuss the properties of a dynamo operat- 
ing in the bulk or in the upper layers of the sun. We en- 
visage a dynamo that operates very much like a classical 
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afl dynamo as it was proposed in the early days of mean 
field dynamo theory. In particular, one may anticipate 
a field that is not strongly fibril, as is indeed confirmed 
by simulations of turbulent dynamos with shear (Bran- 
denburg, Bigazzi, & Subramanian 200f). Furthermore, 
in the uppermost 3 Mm of the sun turbulent downward 
transport is far too strong to let any significant field ap- 
pear at the surface, except in active regions which emerge 
as the result of strong flux segregation into strongly and 
weakly magnetized regions, as demonstrated by large 
aspect ratio magneto-convection simulation (Tao et al. 
1998). 

There are however several new problems. Most impor- 
tant is the fact that the near surface shear layer is rather 
thin, so we may have a problem of too many toroidal 
flux belts if the dynamo was solely confined to the sur- 
face layer. There are other possible problems where only 
a preliminary discussion can be offered. This includes 
potentially observable effects of near surface magnetic ac- 
tivity on the supergranulation. It is conceivable that the 
supergranulation may show significant alignment with 
the mean field. So far, no such anisotropy has been re- 
ported, although we do know that the cell size of the 
normal surface granulation does change with the cycle 
(Houdek et al. 2001). 

The other aspect concerns Joy's law which has success- 
fully been reproduced within the framework of the thin 
flux tube approximation. In the context of a distributed 
dynamo the inclination of bipolar regions is primarily 
controlled by the sense of the latitudinal shear (patches 
at higher latitudes lack behind those at lower latitudes). 
There is also a contribution of a to the tilt (positive 
a produces positively helical fields whose interceptions 
with the surface yield a solar tilt). The latter effect is 
however subdominant, as is seen from a turbulence sim- 
ulation (see the next section). 

We should mention the phenomenon of active zones. 
They constitute patches of recurrent magnetic activ- 
ity over months and sometimes years. Recent inves- 
tigations by Benevolenskaya et al. (1999) showed that 
these patches have different angular velocity at different 
depths, corresponding to the local angular velocity at 
radii between 0.95 R and 0.98 R, suggesting again that 
these magnetic activity complexes are anchored within 
the near-surface shear layer; see Fig.^ an d the gray dots 
in Fig. [2] One may picture these activity complexes as 
more strongly magnetized regions which can only decay 
slowly, possibly because of magnetic helicity conserva- 
tion of perhaps because magnetic fields can have a ten- 
dency to segregate into strongly and weakly magnetized 
regions, as is found in large aspect ratio magnetoconvec- 
tion with imposed field (Tao et al. 1998). 

An important consideration for tachocline dynamos is 
whether the observed emergent flux of 10 24 Mx over the 
full solar cycle can be produced (Galloway & Weiss 1981). 
In the present context we are thinking of mean toroidal 
fields of the order of 300 G, which is about one tenth of 
the strength of the mean field usually envisaged for the 
overshoot layer. However, because of the larger cross- 
sectional surface of, say, (200 Mm) 2 , this will still produce 
the required 10 24 Mx. 

Global simulations of the solar dynamo are becoming 
more advanced. The models of Brun et al. (2004) show 
distributed dynamo action of the type envisaged in this 



paper. However, these models have some characteristic 
properties that are different from the observed solar field. 
Most important is perhaps the relatively large ratio of 
poloidal to toroidal field, which suggests that the effect 
of the differential rotation is not sufficiently prominent 
or, conversely, the effect of the small scale turbulent field 
is too prominent. In the following we discuss a possi- 
ble cause of this and address the question how this may 
change with increasing resolution and larger fluid and 
magnetic Reynolds numbers. 

There are two distinct properties of turbulent dynamo 
action that seem to depend on microscopic viscosity and 
diffusivity. First, the magnetic field in the simulations 
may be dominated by small scale dynamo action (where 
helicity and shear are unimportant). At small mag- 
netic Prandtl number the small scale dynamo is much 
harder to excite (Schekochihin et al. 2004a, Boldyrev 
& Cattaneo 2004, Haugen et al. 2004) and may be- 
come subdominant, allowing the large scale dynamo ef- 
fect to become more prominent. Second, at high mag- 
netic Reynolds numbers the large scale dynamo time 
scale tends to be constrained by magnetic helicity con- 
servation (Brandenburg 2001). Small scale magnetic he- 
licity fluxes throughout the domain become important to 
allow - and even facilitate - large scale dynamo action 
(Blackman & Field 2000, Kleeorin et al. 2000, 2003, Vish- 
niac & Cho 2001, Subramanian & Brandenburg 2004, 
Brandenburg & Sandin 2004). The convective dynamos 
in simulations of Brun et al. (2004), for example, gener- 
ate only a weak mean field. This raises the possibility 
that such dynamos are of a different type than the large 
scale dynamo that is believed to operate in the sun. For 
this reason we now inspect a somewhat different simula- 
tion that lacks convection and stratification, but where 
there is sufficient shear to generate a prominent large 
scale field. We use this simulation to find some guidance 
regarding the question of how fibril the field of the sun 
is and whether it might provide the right conditions for 
bipolar regions to form. 

4. GUIDANCE FROM SIMULATIONS 

The existence of fibril fields (Parker 1982) is crucial in 
a scenario where strong flux tubes rise through the con- 
vection zone to form sunspot pairs at the surface. The 
fibril nature of the field is also the main reason why mag- 
netic buoyancy may be so important. Fibril fields have 
indeed been seen in simulations of forced and convective 
turbulence where mostly a small scale dynamo is in oper- 
ation (Nordlund et al. 1992, Politano, Pouquet, & Sulem 
1995, Brandenburg et al. 1996). However, this picture 
changes in the presence of strong shear, for example in 
the case of accretion disc turbulence (Brandenburg et al. 
1995) or in the case of imposed shear (Brandenburg et 
al. 2001). We begin with a brief description of the model 
and then discuss whether the field is fibril and whether 
it is able to form bipolar regions. 

4.1. Description of the model 

The model discussed in this paper is basically equiv- 
alent to the model studied recently by Brandenburg & 
Sandin (2004), except that no external field is imposed. 
In this model the turbulence is driven by a forcing func- 
tion / that consists of eigenfunctions of the curl operator 
(with wavenumbers 4.5 < kf < 5.5) and of a large scale 
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component with wavenumber ki = 1. The domain is of 
size 7j7r x 2tt x ^ir, representing a cartesian approxima- 
tion (x,y,z) to a sector in the sun between 0° and 30° 
latitude. Thus, (x,y,z) corresponds to (r, (f>, — 6), where 
<t> is longitude and 9 colatitude. The forcing function is 
arranged such that a mean flow of the form 

U = U coskixcoskiz, (1) 

is driven in the meridional plane —tt/2 < k\X < and 
< k\z < 7r/2. In the following we adopt units where 
ki = 1. The equator is assumed to be at z = and the 
outer surface at x = 0. The bottom of the convection 
zone is at x = —tt/2 and z = ir/2 corresponds to the 
latitude of around 30°, i.e. where the surface angular ve- 
locity equals the value in the radiative interior. In the 
plots below we always display nondimensional combina- 
tions by scaling length with k\ and time with u rms fci. 

In this model there is radial shear near the "bottom" of 
what represents the convection zone and latitudinal shear 
in the upper parts. At the level of simplification neces- 
sary to isolate fundamentally new effects, such as the 
shear-current effect of Rogachevskii & Kleeorin (2003, 
2004) we have refrained from modeling the near-surface 
shear layer. Furthermore, curvature effects are ignored 
and no Coriolis force is included, so the 17 x J effect 
of Radler (1969) is absent. With nonhclical forcing, the 
W x J (where W = V x U is the mean vorticity) is how- 
ever a possible effect driving large scale dynamo action. 
Nevertheless, we focus on the morphology of the field in 
the case where the helicity of the forcing is finite and 
negative - consistent with the conditions in the northern 
hemisphere of the sun. 

We assume an isothermal equation of state with sound 
speed c s = const, and solve the continuity, momentum, 
and induction equations in the form 



^ = -cl V lnp + i±* +f + F visc , (3) 
Dt p 

BA 

— = UxB- mo J, (4) 

where p is density, U is velocity, J = V x B/^q is cur- 
rent density, B = V x A is magnetic field expressed in 
terms of the magnetic vector potential, r\ is the magnetic 
diffusivity, and 

F visc = v(V 2 U +±VV -U + 2S- Vlnp) (5) 

is the viscous force where v = const is the kinematic 
viscosity, and Sy = \ (Ui t j + Uj t i) — \8ijUu,k the traceless 
rate of strain tensor. 

We step the equations forward in time by using the 
Pencil Code 1 , which is a high order finite difference 
code (sixth order in space and third order in time) for 
solving the compressible hydromagnetic equations on dis- 
tributed memory machines using the Message Passing 
Interface libraries. The numerical resolution is 128 x 
512 x 128 meshpoints. The boundary conditions are pe- 
riodic in the y direction and stress-free in the x and z 
directions. For the magnetic field we assume perfect con- 
ductor boundary conditions at what corresponds to the 
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Fig. 3. — Evolution of the energies of the total field (B 2 ) 
and of the mean field (B ), in units of B 2 q , for runs with non- 
helical forcing and open or closed boundaries; see the solid and 
dotted lines, respectively. The inset shows a comparison of the 

ratio (B )/(B 2 ) for nonhelical (a = 0) and helical (a > 0) runs. 
For the nonhelical case the run with closed boundaries is also shown 
(dotted line near (B ) / (B 2 ) « 0.07). Note that saturation of the 
large scale field occurs on a dynamical time scale; the resistive time 
scale is given on the upper abscissa. 



base of the convection zone {k\x = —tt/2) and at mid- 
latitudes [k\z = 7r/2). At the equator (z — 0) and at the 
outer surface (x = 0) we assume the magnetic field to be 
normal to the boundaries, i.e. B x n = 0. We refer to 
these boundaries as "open" , because they permit mag- 
netic and current helicity fluxes, as opposed to "closed" 
or perfectly conducting boundaries where B ■ h = with 
no helicity fluxes. 

The strength of the forcing is chosen such that the typ- 
ical rms velocity of the turbulence (without the system- 
atic shear flow) is subsonic with u lnls /c s = 0.1. Viscosity 
and magnetic diffusivity are chosen such that vj-r] = 1 
and R m = u rm s/(?7fcf) ~ 80. In the simulations with neg- 
ative helical forcing, the turbulent part of the velocity 
field is nearly fully helical, i.e. (a; • u)/(kfu^ ms ) « — 1. 
The mean flow is about 5 times stronger than the turbu- 
lent flow, i.e. Uo/u Tms « 5. 

4.2. Growth of the large scale field 

A series of experiments has been performed: with he- 
lical forcing of negative helicity (representative of the 
northern hemisphere; denoted by "a > 0" since the re- 
sulting electromotive force would lead to a positive a 
effect), positive helicity (mainly for comparison, but it 
could be representative of the bottom of the convection 
zone; denoted by "a < 0"), and without helicity (de- 
noted by "a = 0"). In all cases the kinematic growth 
rate is about the same (dlni? rms /dt = 0.02 it rms fcf ). 

In Brandenburg & Sandin (2004) the effect of bound- 
aries was already found to be important: when a per- 
fect conductor condition was used at the equator and 
at the outer surface the resulting a effect was found to 
be suppressed by a factor of ~ 30. In the present case 
we find that near saturation the large scale field remains 
well below equipartition (see the dotted line in Fig. 01 ■ 
With open boundary conditions, near-equipartition field 

strengths can be achieved (B IB 2 w 0.8). Here we de- 



fine B as an average over the y direction (toroidal aver- 
age). Volume averages are denoted by angular brackets. 

In the presence of finite helicity the result is not greatly 
affected; see the inset of Fig. [3] where we show that the 

ratio (B )/(B 2 ) either varies around 0.5 (for a > or 
a < 0), or that it stays around 0.7 (for a = 0). The case 
a < is not shown here, but we refer to Brandenburg et 
al. (2005) for a description of those results. 

So far we have not seen reversals of the field. We 
note, however, that in principle cycles are possible in this 
type of geometry and have indeed be found both in the 
corresponding mean field model (Brandenburg & Sandin 
2004). In the present case the lack of cycles could be con- 
nected with the resulting mean flow that was neglected 
in the mean field calculations. [On the average, however, 
the mean poloidal flow (poleward at the surface) is less 
than 2% of the total mean flow. By comparison, the 
mean poloidal field is about 25% of the total mean field; 
see flOl ] 

Our main conclusion from these simulations is that 
large scale dynamo action can produce equipartition field 
strengths on a dynamical time scale, provided the bound- 
aries are open. The significance of open boundaries is 
that magnetic and current helicity can leave the domain, 
thus preventing the excessive build-up of small scale mag- 
netic helicity before the large scale field has saturated. 

4.3. How fibril is the field? 

Virtually all dynamo simulations (both small scale and 
large scale) show that once the dynamo saturates, the 
typical length scale of the field increases, as measured 
for example by the magnetic Taylor microscale Am 7 where 
a m = 5/io( j2 }/( s2 >; see Schekochihin et al. (2004b) for 
the case of a forced small scale dynamo and Branden- 
burg et al. (1996) for the case of a small scale dynamo 
in convective turbulence. In practice this means that the 
typical scale of the flux structures increases during the 
saturation. However, the orientation of the field is oth- 
erwise still random. The presence of shear together with 
turbulent diffusion has a strong tendency to order the 
field such that it points everywhere in the same direction. 
This tendency has been studied earlier in connection with 
a completely random (incoherent) a effect (Vishniac & 
Brandenburg 1997). This tendency is seen in the present 
simulations as well; see Fig.0J where we show meridional 
cross-sections of the field during kinematic and saturated 
phases. From the simulations presented here we cannot 
support the assumption that in the solar convection zone 
the field will be highly fibril. 

Compared to the sun, there is of course the difference 
that the turbulence is not driven by a body force, but by 
convection. However, it is not clear that this will make 
an important difference. It is also possible that at larger 
magnetic Reynolds numbers there will be a stronger ten- 
dency to produce intense fibrils. However, for the present 
simulations with 128 x 512 x 128 meshpoints the viscos- 
ity is already as small as possible. In fact, the mag- 
netic Reynolds number based on the mesh spacing, Sx, 
of R m kfSx ~ 5, which is a typical value that should not 
be exceeded in these type of simulations. 

Once we abandon the rising flux tube picture, we have 
to think of other ways to produce bipolar regions with 
the right tilt angle. This will be discussed in the next 
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Fig. 4. — Snapshots of the magnetic field B in the meridional 
plane during the kinematic stage (t = fOO turnover times) and 
the saturated stage (t = 1000 turnover times). Vectors in the 
meridional plane are superimposed on a gray scale representation 
of the azimuthal field. The gray scale is symmetric about mid- 
gray shades, so the absence of dark shades (e.g. in the lower panel) 
indicates the absence of negative values. Note the development 
of larger scale structures during the saturated stage with basically 
unidirectional toroidal field. 



section. 

4.4. Formation of bipolar regions 

In the present simulations, because of shear, most of 
the field is in the streamwise direction. On the open 
boundaries, on the other hand, the field can only be 
normal to the boundary and hence B y = 0. However, 
because elsewhere in the interior, the field is mostly 
toroidal, places with significant normal field (B ■ n / 0) 
will be rare. This is also what is seen in the simulations; 
see Fig. were we show "magnetograms" of the normal 
field on the outer surface, x — 0. A meridional cross- 
section of the azimuthally and temporally averaged field, 
B, is shown in Fig.0 Here the identification with a sec- 
tor in a meridional plane of the sun is annotated on the 
axes. Except near the open boundaries, where B y = 0, 
the mean field is mostly into the plane (B y > 0) and is 
accompanied by a right-handed swirl so that B z < on 
the outer surface, x = 0. 

What the magnetograms in Fig. [S] show is basically a 
gray background (corresponding to zero field) with only a 
few patches, some of which come in pairs. Often the pairs 
are connected by a faint "bridge" . In all cases the bipolar 
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Fig. 5. — Gray scale representation of the azimuthally and 
time averaged mean azimuthal field B(x, z) together with vectors 
in the meridional plane. The mean toroidal velocity is shown as 
white contours. The projected positions on the sun are labeled on 
the corresponding axes. Note the equatorward orientation of the 
poloidal field near the outer surface (x = 0). As in Fig. 0] the gray 
scale is symmetric about mid-gray shades, so the absence of dark 
shades indicates that B y > 0. In this run the kinetic helicity is 
negative, corresponding to a > 0. The root mean square of the 
mean poloidal field is ~ 25% of the total mean field. 



regions as well as the bridges are inclined relative to the 
toroidal direction by an inclination angle that is primar- 
ily determined by the latitudinal shear. Any cross-stream 
(i.e. latitudinally oriented) field becomes sheared out and 
intensifies until the structure disappears. The polarity 
depends on the sign of the latitudinal (z) component of 
the field beneath the surface. The phases of maximum 
intensity correspond to times when structures are most 
prominent; see Fig. El Control simulations with opposite 
sign of helicity confirm that the inclination is qualita- 
tively unchanged and that the poloidal field determines 
the orientation of the polarities, not the toroidal field. 
If the tilt was entirely determined by the negative a, it 
would have produced tilt in the opposite sense. In the 
present case of negative kinetic helicity in the northern 
hemisphere (a > 0), as in the upper parts of the solar 
convection zone, the a effect would produce tilt of the 
same sign as the shear. However, as we have seen above, 
the effect of a on the tilt is subdominant in the present 
simulations, where the effect of shear is strong. 

The typical separation of the different polarities corre- 
sponds to about 1-2 eddy scales (= 2'w/kf). Comparing 
to the sun, the pressure scale height at r = 0.95i? is 
12 Mm, so the mixing length and hence the eddy scale is 
~ 20 Mm. Thus, the typical size of the active region in 
the model is ~ 30 Mm. 

We conclude from this section that the tilt of bipolar 
regions depends mainly on the latitudinal differential ro- 
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Fig. 6. — Magnetograms of the radial field B x (y, z, t) at the outer 
surface (x = 0) on the northern hemisphere at times t/r = 480, 
750, and 990, where r = (« r ms%) _1 is the turnover time, and kf is 
the wavenumber corresponding to the energy-carrying scale of the 
turbulence. Light shades correspond to field vectors pointing out 
of the domain, and dark shades correspond to vectors pointing into 
the domain. The elongated rings highlight the positions of bipolar 
regions. Note the clockwise tilt relative to the y (or toroidal) direc- 
tion, and the systematic sequence of polarities (white left and dark 
right) corresponding to B v > 0, which is consistent with Fig. 151 



tation, and that the orientation of the polarities depends 
on the orientation of the latitudinal component of the 
field rather than its azimuthal component. The sense of 
the tilt is thus independent of the sign of a. The simu- 
lations suggest that more or less isolated bipolar regions 
can emerge in a way that is at least as plausible as the 
picture of strong tilted flux tube poking through the sur- 
face from deep underneath. 

5. CONCLUSIONS 

It has long been known that sunspots and other mag- 
netic tracers rotate faster than the photospheric plasma. 
Originally, this was taken as evidence that the sun must 
rotate faster in the interior (Golub et al. 1981). In fact, it 
was believed that the angular velocity of sunspots agrees 
with the local angular velocity at the depth where the 
sunspots are anchored. Since the mid-eighties this idea 
became largely discarded on the grounds that helioseis- 
mology began to show angular velocity contours that are 
nearly spoke-like and that the only location of radial 
shear was the bottom of the convection zone. From a 
dynamo theorist's point of view this result, together with 
the already popular idea that the solar dynamo should 
operate at or below the bottom of the convection zone 
(Table HJ, meant that at least the issue of the location of 
the dynamo had been settled. This picture was in princi- 
ple quite appealing and it became particularly attractive 
in combination with the subsequent finding that the right 
tilt angles can be obtained when the fields at the bottom 
of the convection zone are of the order of 10 5 G (D'Silva 
& Choudhuri 1993). 

Two important results have emerged since then. First, 
downward pumping tends to be a strong effect that 
can overcome magnetic buoyancy up to fairly large field 
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strengths. Second, helioseismology has now revealed the 
presence of a near-surface shear layer that is stronger 
and more prominent than indicated by the early results 
of helioseismology. In the tachocline, by comparison, the 
radial shear layer is rather weak at 30° latitude (which is 
where sunspots emerge in the beginning of the cycle) and 
extremely strong at the poles (where the magnetic activ- 
ity is weak) . It appears that the impact of these findings 
on the solar dynamo paradigm ought to be reconsidered. 
The purpose of the present paper was therefore to present 
the arguments for and against tachocline dynamos versus 
distributed dynamos that are possibly strongly affected 
or shaped by the near-surface shear layer. The reason 
we use word "shaped" is to indicate that the dynamo is 
likely to operate in the entire convection zone, and not 
only in the near-surface shear layer. 

The idea that sunspots might be anchored at a depth 
of 0.95 solar radii has however a problem. The helioseis- 
mologically determined angular velocity at that depth 
is still ~ 4 nHz slower than that of the very youngest 
sunspots. This corresponds to a velocity difference of 
~ 20m/s. In principle, if the profile of such enhanced 
angular velocity is sufficiently localized, one might argue 
that the spatial resolution of helioseismology was still in- 
sufficient. Alternatively, and perhaps more likely, there 
may still be some not yet understood mechanism causing 
newly emerging flux to rotate slightly faster. 

An important part of the flux tube paradigm is the 
idea that strong flux tubes emerge from deep underneath 
and form bipolar regions and sunspot pairs as they reach 
the surface. This picture relies entirely on the thin flux 
tube approximation, which may have its own difficulties 
(Dorch & Nordlund 1998, Wissink et al. 2000). However, 
as demonstrated in 21 and through Fig. EI a distributed 
dynamo is quite able to produce bipolar regions with 
plausible tilt angles. As discussed at the end of M4.4I 
the orientation of the tilt is controlled by the latitudinal 
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shear. The orientation of the polarities is determined by 
the direction of the poloidal field, and not the azimuthal 
field. 

Since this model lacks convection and stratification, 
both turbulence and shear have to be produced by body 
forces. Nevertheless, the model is fully self-consistent 
and not subject to approximations, such as the thin 
flux tube approximation. Obviously, an important next 
step should be to include convection and stratification. 
Equally important is the implementation of a more re- 
alistic outer boundary condition, possibly allowing for 
the development of coronal mass ejections that might be 
necessary for carrying small scale magnetic and current 
helicities away from the dynamo. Proper modeling of 
coronal mass ejections might require the use of spherical 
geometry. A lack of magnetic and current helicity fluxes 
out of the domain would prevent the dynamo from op- 
erating on a dynamical time scale (Blackman & Field 
2000). On the other hand, some degree of throttling of 
the helicity flux might actually occur in the sun. This 
might explain why the solar cycle period tends to be 
about 10 times longer than what is suggested by stan- 
dard mean field models (Kohler 1973) using canonical es- 
timates for the turbulent diffusivity (Krivodubskii 1984). 
Obviously, a reasonably accurate theory of helicity fluxes 
is required before this question can be addressed in mean 
field calculations. 
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